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Abstract 

The heavy fermion system CeCu6-xAg x is studied at its antiferromagnetic quantum critical 
point, x c = 0.2, by low temperature (T > 50 mK) specific heat, C(T), and volume thermal expan- 
sion, (3(T), measurements. Whereas C/T oc log(To/T) would be compatible with the predictions 
of the itinerant spin-density-wave (SDW) theory for two-dimensional critical spin fluctuations, 
(3{T)/T and the Griineisen ratio, T(T) oc (5/C, diverge much weaker than expected, in strong con- 
trast to this model. Both C and (3, plotted as a function of the reduced temperature t = T/Tq with 
To = 4.6 K are similar to what was observed for YbRli2(Sio.95Geo.o5)2 (To = 23.3 K) indicating a 
striking discrepancy to the SDW prediction in both systems. 

PACS numbers: 71.10.HF,71.27.+a 



The nature of the quantum critical point (QCP) at which long-range magnetic order 
gradually develops in heavy-fermion (HF) systems has been the focus of research activities 
during past years Q]. In paramagnetic systems like CeCug, a Fermi surface formed by heavy 
quasiparticles has been detected, indicating an itinerant 4/-electron state Q|. The central 
question is how these heavy quasiparticles evolve if these materials are tuned into a long- 
range magnetically ordered state. In the traditional picture , the quasiparticles retain 
their itinerant character and form a spin-density wave (SDW) type of antiferromagnetic (AF) 
ordering. For this Hertz-Millis-Moriya theory, in the following refered as itinerant scenario, 
a mean-field type of quantum critical behavior arises. Recent experiments have shown that 

no 

at least in some HF systems this picture fails p, [7|. Consequently, a new type of QCP 
has been proposed at which the quasiparticles break-up into their components: conduction 
electrons and local 4/-moments forming magnetic order. This locally- critical scenario arises 



due to the destruction of the Kondo resonance at the QCP 



aa. 



In CeCu6, the weakening of the 4 /-conduction electron hybridization, J, induced by a 
lattice expansion in CeCu 6 _ x JVL with M=Au (x c = 0.1) Q, Ag (x c = 0.2) 11, Q, Pd 
(x c = 0.05) and Pt (x c = 0.1) [l^J drives the system at a critical concentration x c towards a 



QCP beyond which long-range AF order appears. Remarkably, the non-Fermi liquid (NFL) 
behavior observed at the respective x c in the different systems does not depend on which 
dopand M has been used ^J. In fact, the specific heat coefficient of all of these different sys- 
tems is very similar (within ~ 10% below 3 K) and shows a characteristic C/T oc log(To/T) 
dependence over nearly two decades in T down to 50 mK. Thus, a common origin of the 
quantum critical state at x = x c for the different systems CeC^-xM^ is likely. The log- 
arithmic divergence of the quasiparticle mass would only be compatible with the itinerant 
scenario if strongly anisotropic, two-dimensional (2D) critical spin fluctuations would lead 
to singular scattering at the whole Fermi surface Q|. Indeed, inelastic neutron scattering 
experiments on CeCus.gAuo.i, revealed rod-like structures of high intensity in g-space trans- 
lating to quasi-2D fluctuations in real space jl^. The 2D SDW picture, however, cannot 
explain that both at and far away from the wave vector of the nearby AF order the neutron 
and magnetization data exhibit E/T and B/T scaling with an anomalous exponent a < 1 
The momentum independence in the critical response led to the proposal of local criti- 
cality 0,0,0]. So far it is not clear, which kind of mass divergence would occur in such a 
scenario. At the AF QCP in YbRli2(Sio.95Geo.o5)2 a stronger than logarithmic mass diver- 
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gence has been observed at very low temperatures, incompatible with the itinerant scenario 
7]. Furthermore, local moments have been detected in the low- frequency bulk susceptibility 
7], as well as in microscopic magnetic measurements 16, Q], indicating local criticality in 



this system. Comparing the title Ce system with this Yb system, two important questions 
arise: i) is there a similar discrepancy in thermodynamic properties with respect to the pre- 
dictions of the itinerant scenario for the Ce system as found 7J for YbRh 2 (Sio.95Geo.os)2 and 
ii) do both materials show similar mass divergences? This Letter addresses these questions. 

Very recently, it has been shown that the thermal expansion, (3 = V~ 1 (dV/ dT) (V: sample 
volume), compared to the specific heat, behavior is much more singular in the approach to 
the QCP. Complementary to the specific heat, being related to the temperature dependence 
of the entropy, the thermal expansion probes the pressure dependence of the entropy and 
;hus should become singular in the approach of a (pressure-driven) QCP. Scaling analysis 
3| revealed that the Gruneisen ratio r oc (3/C can be used as a highly sensitive probe of 
quantum criticality because it has to diverge at any QCP. The critical exponent e in the 
divergence V oc 1/T e is given by e = 1/vz with u, the critical exponent for the correlation 
length, £ oc \r\ v (r: distance from the QCP) and z, the dynamical critical exponent in the 
divergence of the correlation time, r c oc £ z [18]. For a 3D AF QCP ^ the itinerant scenario 
predicts v = 1/2 and z = 2 yielding e = 1. Thus a study of the Gruneisen ratio can serve as 
a basis for a detailed assessment of the validity of the SDW picture in different systems. 

The first-ever observation of a Gruneisen ratio divergence has been made on the two 
HF systems CeNi 2 Ge 2 and YbRh 2 (Sio.95Ge .o5)2, both located at AF QCPs Q- For the 
former system e = 1 has been found within 50 niK < T < 3 K in accordance with the 
itinerant scenario. For YbRh 2 (Si .95Ge .o5)2, on the other hand, the measured Gruneisen 
exponent is fractional: e = 0.7 ± 0.1 for T < 0.8 K j^J. This cannot be explained by the 
itinerant theory, but is qualitatively consistent with the locally quantum critical picture. As 
this scenario has first been proposed on the basis of the neutron scattering experiments on 
CeC\iQ- x M x (M=Au) f|, it would be highly desirable, to determine the critical exponent of 
the Gruneisen ratio divergence in this system at the critical concentration x c . 

Below, we show first that the QCP in CeCu6_ x Ag a; (x c = 0.2) is phenomenologically re- 
lated to that in CeCug-xAu^ (x c = 0.1) and afterwards investigate the nature of this QCP by 
means of a Gruneisen analysis. Polycrystalline samples of CeCug-xAg^ with 0.09 < x < 1.2 
were prepared via arc melting together stoichiometric amounts of high-purity elements under 



3 



purified argon atmosphere |12l.l2lLl22|. They were analyzed by X-ray powder diffraction and 
found to be single phase and with the proper orthorhombic structure. In contrast to the Au 
system where for x < 1 only one of the five different Cu sites is occupied by the dopand, 
the Ag atoms are randomly distributed on all Cu sites 2^|. The linear enhancement of the 
b- and c-lattice parameters is observed to be of similar size as in the Au system. On the 
other hand, at x = 1 the a-parameter enhancement is 4% larger compared to that in the 
Au system. Specific heat and thermal expansion measurements were performed in dilution 



refrigerators with the aid of a thermal relaxation calorimeter 1 121 and capacitive dilatometer 

n 

|20l |. respectively. 

Figure 1 displays the low-temperature specific heat divided by temperature, C(T)/T, of 
various CeCu6-iAg x . samples on a logarithmic temperature scale. Long-range AF order is 
observed for x > 0.3 and is manifested by broadened jumps in C(T)/T. The inset shows 
Tjv(x) as determined by (entropy-conserving) equal areas constructions and the maximum 
of the derivative dp(T)/dT in corresponding electrical resistivity measurements j^j]. Ex- 
trapolation of T/v to zero temperature yields a critical concentration x c = 0.2. This value 
is in accordance with the extrapolation of the entropy at the AF phase transition, S(Tn), 
measuring the size of the ordered moment, towards zero. Thus, the critical concentration is 
twice as large as in the case of the CeC^-zAu^. system although the relative volume expan- 
sion by Ag doping is slightly larger as in the Au case. It is unlikely that the larger value of 
x c is caused by stronger site disorder in the Ag system. Another possibility is that the Ag 
ions adopt a different valence state than the Au ions j^jj. As the volume of Ag 1+ ions is 
larger than that of Au 2+ ions this would explain the relatively larger volume expansion for 
Ag, compared to Au. Furthermore, a change in the electronic configuration might increase 
the density of states at the Fermi level. This would effectively enhance J and stabilize the 
paramagnetic ground state leading to a higher x c for the Ag system |2l[ . 

For the CeCus.sAgo^ sample which is located right at the QCP, the specific heat coefficient 
is very similar to that in CeCus.gAuo.i and, at zero magnetic field, diverges logarithmically 
between 50 mK and 2.5 K (Figure 2a). The application of magnetic fields leads to the 
recovery of Landau Fermi-liquid behavior, C/T = 7(-B), at lowest temperatures as observed 
in many other NFL systems We observe a logarithmic divergence 7(5) oc — logB (not 
shown) which represents a clear signature of a QCP at B = 0. The zero-field temperature 
dependence, as well as the field dependence of 7(-B), assuming a linear relation between the 
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magnetic field and the control parameter, would be compatible with the itinerant theory for 
an AF QCP, if one assumes truely 2D critical spin fluctuations. 

We now turn to the thermal expansion, measured along three perpendicular orientations 
on the same polycrystal studied by specific heat (Figure 2b). The observed anisotropy 
is caused by the texture of the polycrystalline sample. The volume expansion coefficient 
(3 is given by the sum of the three linear expansion coefficients on all showing a simi- 
lar temperature dependence. Upon cooling to the lowest temperatures, f3(T)/T increases 
strongly and diverges logarithmically for T < 0.8 K. Although the observed (3/T diver- 
gence is steeper than in C/T, indicating that the NFL behavior is caused by a QCP, 
it is much weaker than the temperature dependence expected in the 2D-SDW scenario, 
f3(T)/T = a + ai /T + a 2 /T log log (T /T) Q. 

Figure 3 displays the temperature dependence of the Griineisen ratio T(T) calculated from 
the specific heat and thermal expansion data shown in Figure 2. In the entire temperature 
range, the divergence is weaker than T oc 1/T and thus incompatible with the predictions 
of the itinerant scenario for both 3D or 2D critical spin fluctuations [3]. As shown in the 
inset, T(T) roughly follows a logarithmic increase upon cooling from below 1 K. 

In the last part of this paper, we compare the NFL behavior in Ce Cu^ sAg n ? with that 
observed in YbRh 2 (Sio.95Geo.o5 )2 BQ 

As discussed by Sereni et al. [24j for a substantial 



number of HF systems located near the magnetic instability the logarithmic increase of the 
specific heat coefficient C/T involves a similar amount of entropy (about 55% of i?ln2). 
On a reduced temperature scale t = T/T , where the spin- fluctuation temperature T is 
assumed to be close to the Kondo temperature Tk determined, e.g. from the residual 
quasielastic neutron line width, the specific heat coefficient is described by C/t = —D \ogt + 
ET with D = 7.2 Jmol^K" 1 Q Here the parameter E is small or negligible for systems 
showing a low Kondo temperature as well as excited crystal electric field states that are well 
separated from the ground state. In Figure 4, we show that our specific heat data agree 
perfectly with the scaling function using To = 4.6 K (see straight line). Similar scaling 
behavior has been reported for the CeC\iQ- x An x series as well |24j. The specific heat data 
of YbRh 2 (Sio.95Ge .o5)2 collapse for t > 0.01 on the scaling function with T = 23.3 K (see 
Figure 4). Indeed, this value agrees well with the single-ion Kondo temperature Tk ~ 25 K 
for this system [2|. For t < 0.01 the specific heat coefficient of YbRh 2 (Sio.95Ge .o5 2 shows a 
pronounced upturn that is not a precursor of the tiny AF order at T/v = 20 mK 7] . Above 
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50 mK, the zero-field data of both specific heat and thermal expansion are identical to their 
counterparts at the critical field B c = 0.027 T |2(3]. Since the stronger than logarithmic 
mass divergence in this compound is at strong variance to the predictions of the itinerant 
model, it would be highly desirable to probe the specific heat of the Ce-based systems at 
t < 0.01. For CeCu5.gAgo.2, this corresponds to the temperature range below 50 mK which 
has, unfortunately, not been investigated so far. 

Finally, we turn to the scaling behavior observed in the thermal expansion. In contrast 
to the specific heat, the thermal expansion of CeCu 5 . 8 Ag 2 deviates strongly from the pre- 
dictions of the itinerant theory for T > 50 mK. In Figure 5 we demonstrate that, using the 
same T n values obtained from the specific heat scaling, the thermal expansion shows scaling 
behavior for both systems as well. As mentioned before, the SDW theory in the 2D case 
would predict a very strong divergence (3/T oc T _1 loglogT. By contrast, a logarithmic 
increase is observed upon cooling over more than one decade in t for 0.015 < t < 0.25, 
followed by a 1/t divergence in the Yb system. Already the weak logarithmic divergence 
above 0.015£ rules out the 2D SDW scenario for the QCP in both systems. 

In summary, the thermodynamic behavior in CeCu5.8Ago.2, which is located right at an 
AF QCP, has been studied by measurements of the specific heat and the thermal expansion. 
The observed logarithmic temperature dependence of C /T would be compatible with the 
itinerant theory assuming the presence of 2D critical spinfluctuations. Such a scenario, can 
however clearly be excluded from the analysis of the thermal expansion and Griineisen ratio 
data. Both properties diverge towards zero temperature much weaker than predicted by the 
itinerant theory. On the reduced temperature scale T/Tq, with To ~ Tk, the Kondo tempera- 
ture, one finds universal behavior in the thermodynamic properties of both CeCus.gAgn^ and 
YbRh? (Sin QsGen nO ? for which latter material a locally critical QCP had been highlighted 
recently p, |l7|, boj . 
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FIG. 1: Specific heat as C/T vs T (on a logarithmic scale) for different CeC\\§- x A.g x polycrystals. 
Inset shows evolution of antiferromagneticphase transition temperature T/v vs x as derived from 
specific heat (squares: this study, circles [llj]) and electrical resistivity (triangles [2l|]) results. 



8 




0.05 0.1 



T(K) 



T(K) 



FIG. 2: (a) Electronic specific heat of CeCu5.sAgo.2 as C e \/T vs T (on a logarithmic scale) for 
B = and differing magnetic fields. At B > 0, C e i is obtained after subtraction of the Cu nuclear 
specific heat contribution C n oc B 2 /T 2 j^|. (b) Volume thermal expansion coefficient (3 of the 
same sample studied in (a) as 0/T vs logT. (3 = a\ + 02 + «3 with a% being the linear thermal 
expansion coefficients along the three perpendicular directions of the sample. Solid lines indicate 
logarithmic temperature dependences in C e i{T)/T and (5{T)/T. 
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FIG. 3: Temperature dependence of Griineisen-ratio T = V-moil^T • P/C wit 
6.37 • 10~ 5 m 3 mol and isothermal compressibility kt = 1 • 10~ n Pa -1 
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23j. Inset shows same 



data on a logarithmic temperature scale. Solid line represents — log(T) dependence. 
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FIG. 4: Scaling of the low-T specific heat of CeCus.sAgn 2 and YbRh2(Sio.95Geo.o5)2 |2L I2UJ as 
(C/T — E)Tq vs log(T/To) according to Sereni et al. [24|. For CeCu5.sAgo.2, Tq = 4.6 K and 
E = 0.105 J/K 2 mol, for YbRh2(Si .9 5 Geo.o5)2, T Q = 23.3 K and E = 0.066 J/K 2 mol. Solid line 



represents scaling function 7.2 log (Tq/T) observed in several other HF systems 
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FIG. 5: Volume thermal expansion as (5/T vs normalized temperature T /Tn (on a logarithmic 
scale) for CeCu 5 . 8 Ag .2 (T = 4.6 K, right axis) and YbRh 2 (Sio.9 5 Geo.o5)2 [2o| (T = 23.3 K, left 
axis). Parallel solid lines indicate logarithmic temperature dependence of (5{T)/T. 
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